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The entanglement dynamics of two independent qubits each embedded in a structured environ- 
ment under conditions of inhibition of spontaneous emission is analyzed, showing entanglement 
trapping. We demonstrate that entanglement trapping can be used efficiently to prevent entangle- 
ment sudden death. For the case of realistic photonic band-gap materials, we show that high values 
of entanglement trapping can be achieved. This result is of both fundamental and applicative inter- 
est since it provides a physical situation where the entanglement can be preserved and manipulated, 
e.g. by Stark-shifting the qubit transition frequency outside and inside the gap. 
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Entanglement preservation is an important challenge 
in quantum information and computation technologies 
Realistic quantum systems are affected by decoher- 
ence and entanglement losses because of the unavoidable 
interaction with their environments For example in 
Markovian (memoryless) environments, in spite of an ex- 
ponential decay of the single qubit coherence, the en- 
tanglement between two qubits may completely disap- 
pear at a finite time [1, . This phenomenon, known as 
"entanglement sudden death" and proven to occur in a 
quantum optics experiment in turn limits the time 
when entanglement could be exploited for practical pur- 
poses. It is therefore of interest to examine the possibility 
to preserve entanglement. In the case of environments 
with memory (non-Markovian), such as imperfect cavi- 
ties supporting a mode resonant with the atomic transi- 
tion frequency, revivals of two-qubit entanglement have 
been found [a, SSl- These revivals, although effectively 
extending the possible usage time of entanglement, de- 
crease with time and eventually disappear after a certain 
critical time. Moreover, when the qubits interact with a 
common environment, it has been shown that entangle- 
ment can be preserved by means of the quantum Zeno 
effect [il. 

In this paper we continue the investigation on physi- 
cal systems and physical effects that may lead to effec- 
tive long time entanglement protection. Since entangle- 
ment evolution and population decay have been previ- 
ously shown to be related Q, one is led to investigate 
situations where population trapping occurs. This can 
happen in structured environments where the density of 
states presents a dip which can inhibit spontaneous emis- 
sion in the region of the dip Among realistic physi- 
cal situations, this effect is known to occur in photonic 
band-gap (PBG) materials [13] ■ Entanglement can be 
generated when a pair of atoms near-resonantly coupled 
to the edg e of a PBG present direct dipole-dipole inter- 
action [ll|. On the other hand, a way to produce entan- 
gled independent atoms in PBG materials is to consider 



a three-dimensional photonic crystal single-mode cavity 
with a sufficiently high-quality Q factor where Rydberg 
atoms can freely travel through the connected void re- 
gions [l^ . The atoms exchange photons with the cavity, 
represented by a single defect mode of the crystal res- 
onant with the atomic transition frequency, which acts 
as an atomic entanglement catalyst. In this situation, 
the atoms can be considered independent because their 
transition frequency is well inside the PBG and their rel- 
ative distance is large enough to make the dipole-dipole 
interaction negligible. 

The main aim of this work is to analyze if and to what 
extent entanglement initially present may be preserved 
for a general class of states of noninteracting pairs of 
atoms. Finally, we also analyze how such a system may 
be realistically manipulated to perform quantum compu- 
tation and quantum information processing. 

Our system consists of an entangled couple of iden- 
tical independent two-level atoms (qubits), with states 
|0),|1) and transition frequency Wq, each embedded in 
a PBG material at zero temperature. The atoms are 
assumed to be sufficiently separated in order to neglect 
the dipole-dipole interaction. Under this condition it has 
been shown that the two-qubit dynamics can be obtained 
by the knowledge of that of the single qubit Each 
atom is coupled to the surrounding radiation field in a 



three-dimensional periodic dielectric [13[. In this case, 
the Hamiltonian of the single atom-environment system 
can be put in the form 



ujk - uJo)blbk + igk {o--bl - a+bk 



(1) 



where loq is the transition frequency and a± are the rais- 
ing and lowering operators of the qubit, bk are the 
creation and annihilation operators and gk the coupling 
constant of the mode k with frequency ujk 
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Here d and are the absolute value and unit vector 
of the atomic dipole moment, V is the sample volume, 
e^, = Gfe j are the two transverse (polarization) unit vec- 
tors, and Co the Coulomb constant. The Hamiltonian of 
Eq. dll) is physically equivalent to the standard Hamil- 
tonian of a two-level atom interacting with the radiation 
field and it is obtainable from the last one by unitary 
transformations 
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We consider here a model of electromagnetic waves 
in a three-dimensional periodic dielectric where the pho- 
ton dispersion relation ujk is chosen to be isotropic. By 
symmetrizing tok, one produces PBGs at the spheres 
|k| = tott/L (m = 1, 2, . . .), where L is the lattice con- 
stant. In such ideal photonic crystals, a PEG is the 
frequency range over which the local density of electro- 
magnetic states and the decay rate of the atomic popu- 
lation of the excited state vanish. Near the band-gap 
edges the density of states becomes singular [13|, the 



atom-field interaction becomes strong and one can ex- 
pect modifications to the spontaneous emission decay. In 
terms of the detuning S ^ ujq — lUc between the atomic 
transition frequency ujq and the band edge frequency Wc, 
the atomic population of the excited state is given by 
P{t) =_P(0)|(?(t)p where the quantity q{t) may be writ- 
ten as [l4| 



1 - $ 



q{t) = 2aixie'^"i*+'*'* + ai{x2 + ^xl)e^^-*+'^' 
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cf / Gneohc^ , ^{x) is the error function 



1,2,3) with 

xi = (A++A_)e*("/4\ 

X3 = (A+p5^/6)+A_e-'('^/6))e^(3V4)^ 



and A± = [1/2 ± ^/TTIP/2W]^^^ . In the case when 
the environment is at zero temperature and the atom is 
initially in a general superposition state of its two lev- 
els, the single-qubit density matrix evolution can be ex- 
pressed as y 



pfo(0)g(i) 



p'oMq*{t) pfo(0) + pfi(0)(l-|<z(t)|^ 



(4) 

The expression for q{t) of Eq. ([3]) thus allows to obtain 
the single qubit dynamics according to Eq. (|4]) and then 
to determine the two-qubit dynamics by exploiting a pro- 
cedure that relies on the knowledge of the single-qubit 
reduced dynamics Q. In particular, in the standard ba- 
sis B = {|1) ^ |11),|2) ^ |10),|3) ^ |01),|4) ^ |00)}, 



the diagonal elements of the two-qubit reduced density 
matrix p{t) at time t are 

Pii(i) = pimmw 

P22(t) - Pll(0)k(i)P(l-k(i)n+P22(0)k(t)P, 

- Pii(o)|<zWP(i-kWn + P33(o)k(t)p, 



Pii{t) = I - [pilit) + p22{t) + P33{t)], 



(5) 



and the non-diagonal elements 
Pi2{t) = pi2{0)q{tMt)\\ Pi3(i)-Pi3(0)gWk(i)l', 

Puit) = P14(0)(Z(^)^ P23W =P23(0)|(z(i)l', 
P24{t) = pi3(0)q(t)(l- |g(i)|')+p24(0)g(0, 
P34W = Pi2(0)<z(t)(l - + P34(0)gW, (6) 

with Pij{t) — Pji{t), p{t) being a hermitian matrix. Note 
that Eqs. ([5|) and ([6|) permit to obtain the two-qubit den- 
sity matrix evolution for any initial state. 

To analyze the two- qub it entanglement dynamics we 
use the concurrence C [16|, which attains its maximum 
value 1 for maximally entangled states and vanishes for 
separable states. We consider as initial states the Bell- 
like states 

1$) =a|01) -|-e*''/3|10), |*) = a|00) + e*''/3|ll), (7) 

where a, (3 are real and + 0^ = 1. For these initial 
conditions the concurrence at time t is given respectively 
by 



C^{t) = 2max{0,av/l -a^kWP}, 
C^{t) = 2max{0, -a^kWP 

x[a- v/l-a2(l-k(i)|2)]}. (8) 

For initial Bell (maximally entangled) states, |$) = 
(|01) ± |10)) /^/2 and 1^-) = (|00) ± |11)) /\/2, the con- 
currences at the time t take the simple form 



C^{t) = \q{t)\\ C^{t) = \q{t)\' 



(9) 



independently on the relative phase 7. The previous 
equation is relevant because it shows the direct link be- 
tween the behavior of excited state population of single 
atom P{t) and that of concurrence. 

The concurrences C^{t) and C^{t) of Eq. [HI corre- 
sponding to the initial Bell states |$) and I^E*) respec- 
tively, are shown in Fig. [T] and Fig. [5] as a function of the 
dimensionless time (3t and for various values of detuning 
5, 5 < Q corresponding to the case when the two-level 
atom frequency ojq is inside the band gap. The figures 
display that entanglement trapping is possible. The con- 
currence is nearer to its maximum value 1 for atomic 
frequencies at larger distances from the band edge and 
deeper inside the gap. From Fig. [T] and Fig. [5] no quali- 
tative difference appears between the evolutions of con- 
currences for the two initial states I*!") and l^*). A small 
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FIG. 1: Concurrence C#(t) as a function of /3t starting from the 
initial Bell-like state I*} = (|01> ± |10>) /\/2 for various values of 
detuning from photonic band edge: 5 = —10/3 (solid curve), S = 
—4/3 (dotted curve), 5 = —f) (long-short-dashed curve), 5 = (long- 
short-short-dashed curve), S = l3 (long-dashed curve), S = 10/3 
(short-dashed curve). 
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FIG. 3: Asymptotic values of concurrences C^{t = oo) (solid line) 
and Cqi{t = oo) (dashed line) as a function of 5/13, starting re- 
spectively from the initial Bell state |<I>) = (|01) ± |10)) /V2 and 
|*> = (|00>±|11»/V2. 
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FIG. 2: Concurrence Cisi{t) as a function of (3t starting from the 
initial Bell-like state I*} = (|00) ± |11>) /V2 for various values of 
detuning from photonic band edge: 5 = —10/3 (solid curve), 5 = 
—4/3 (dotted curve), 5 = —f) (long-short-dashed curve), 5 = (long- 
short-short-dashed curve), S = f} (long-dashed curve), S = 10/3 
(short-dashed curve). 



quantitative difference can be appreciated in the asymp- 
totic values of concurrences C$((X)) and C^ii{oo), with 
C,i,(oo) being shghtly lower than C$((X)). The asymp- 
totic values of concurrences for t — s- oo can be obtained by 
the steady-state atomic population in the excited state 
|(7sp — limi^oo k(OP = 4|aia;ip 13] on the basis of 
Eq. ([9]). In Fig. [3] C<i.(oo) and C*(oo) are plotted as a 
function of the ratio (5//3. We see that, when 5/ [3 < 0, 
these asymptotic values differ from zero and decrease 
very rapidly moving near the edge. However, also for 
a small range of (5//3 > 0, they are not zero. This can be 
linked to the presence of a photon-atom bound dressed 
state in the single qubit dynamics . It is evident that 
the more the atomic transition frequency is far from the 
band edge inside the gap < 0), the higher is the 

preserved entanglement. This figure confirms that the 
preserved entanglement starting from the Bell state |\[') 



is smaller than that relative to the initial state |$). As 
said, entanglement trapping is achievable even when loq 
lies outside the band gap near the edge, although with 
small asymptotic values (< 0.4). Exploiting Eq. ([8]) it is 
possible to study the behavior of concurrence in terms 
of the initial degree of entanglement (represented by a). 
The result is that, as one might expect, smaller initial 
entanglement leads to lower asymptotic value of concur- 
rence. Entanglement sudden death and its consequent 
limit on the entanglement usage can be thus prevented 
in such systems. 

The results here obtained are valid for ideal PBG ma- 
terials, thus one may ask to what extent they hold for 
real crystals. In fact, in real crystals with finite dimen- 
sions a pseudogap is typically obtained where the density 
of states is strongly smaller than that of free space but 
is not exactly zero. In such gaps, however, the lowering 
of the local density of states can be so large that the 
spontaneous decay of an excited emitter is effectively in- 
hibited so that memory and coherent control effects are 
admitted and similar to those occurring in an ideal PBG 
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igj l . This would correspond, in view of the rela- 
tion between population and entanglement expressed by 
Eqs. (IS]) and to inhibition of entanglement decay. 

Another aspect that must be considered is the possi- 
bility of performing, in the conditions here considered, 
two-qubit operations for quantum computing. When the 
qubit transition frequencies are inside the gap, the ef- 
fect of any resonant external signal is impeded. Thus, 
a way to bring the qubit frequencies outside the gap is 
required for the time necessary to externally manipulate 
the qubits. This could be realized, e.g., by means of Stark 
shifting the frequencies with a static electric field. The 
duration of this external control should be short enough 
to limit the decay of entanglement. For example, suppose 
that two Rb Rydberg atoms, having transition frequency 
luq ^ 50 GHz [iOl, are embedded in a photonic crys- 
tal built in such a way to have PBG frequencies in the 
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range of GHz 2l|. For such atoms the dipole moment 
is d - 2 X 10"^^ Cm [l3|, that in turn gives /? - 20 
KHz, while the typical Stark shifts are A ~ 200 KHz 
[20]. This implies a dimensionless shift with respect to 
/? of the order A//3 ~ 10. On the basis of the results 
displayed in Fig. O this shift is therefore large enough to 
move the atomic transition frequency, for example, from 
a ratio 5/(3 = —5 inside the gap, corresponding to an 
entanglement trapping value ~ 0.9, to a ratio 6/(3 = 5 
outside the gap, where the atoms can be manipulated by 
an external signal. 

In many experiments the qubits are mimicked by "ar- 
tificial atoms" consisting in quantum dot structures em- 
bedded in the solid fraction of the PEG material 1^, 2^ ■ 
Quantum dots permit the coherent manipulation of a 
single localized quantum system with the technological 
advantages of solid-state systems. In fact, the coherent 
control of an exciton wave function in a quantum dot, 
namely the manipulation of the relative phases of the 
eigenstates in a quantum superposition, is experimentally 
achievable [23|. Also in the quantum dots scenario one 
can tune the energy levels by means of the Stark effect 
with typical shifts 
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the qubit transition frequencies outside the gap for a time 
necessary to the desired two-qubit operation. We have 
pointed out that Rydberg atoms and quantum dots in 
photonic crystals may be suitable to this purpose. These 
considerations clearly highlight the potential of struc- 
tured environments as a promising background in the 
development of new quantum technologies. 
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quired conditions for qubit coherent manipulations could 
be accomplished. These features make quantum dots in- 
corporated in a PEG material good candidates for the 
implementation of various schemes for quantum compu- 
tation and coherent information processing. 

In the context of Rydberg atoms traveling through 
the connected void regions of a photonic crystal, oppor- 
tune atom-cavity interactions allow to perform quantum 
logic gates, as dual-rail Hadamard and NOT gates [l3 |. 
and also other quantum information procedures based 
on cavity QED seem to be achievable. After the gate 
operations, when the atoms are outside the cavity in the 
periodic region of the crystal (PEG), the entanglement 
between the two qubits can be preserved for a very long 
time due to the entanglement trapping phenomenon we 
have demonstrated in this paper. 

In conclusion, in this paper we have presented a new 
phenomenon, namely entanglement trapping in struc- 
tured reservoirs, and we have shown that it can be used 
to prevent entanglement sudden death. We have investi- 
gated the connection between entanglement trapping and 
the well known phenomenon of population trapping oc- 
curring in photonic band gap structures. We have shown 
that the entanglement, as measured by concurrence, ini- 
tially present between independent qubits embedded in 
structured reservoirs, as PEG materials, can be preserved 
depending on the "position" of the qubit transition fre- 
quencies inside the band gap. The qubits may be coher- 
ently manipulated by an external control field to perform 
some quantum computation or information protocol. A 
possible way to achieve this is by a Stark shift that moves [25] 



[1 
[2: 

[s: 
[4: 

[s: 
[6: 

[7] 

[8 

[9 
[10 

[11 

[12 

[13 
[14 
[15 



[16: 

[17 

[is: 

[19 

[20 

[21 
[22 
[23 
[24 



Electronic 



address: 



|lofranco@fisica. unipa.it [ 



URL: http : //www.f isica.unipa. it/~lof ranco 

M. A. Nielsen and I. L. Cliuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
2000). 

H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, Oxford, 
New York, 2002). 

L. Diosi, Lect. Notes Phys. 622, 157 (2003). 

T. Yu and J. H. Eberly, Phys. Rev. Lett. 93, 140404 

(2004). 

M. P. Almeida et al.. Science 316, 579 (2007). 

B. Bellomo, R. Lo Franco, and G. Compagno, Phys. Rev. 

Lett. 99, 160502 (2007). 

B. Bellomo, R. Lo Franco, and G. Compagno, Phys. Rev. 
A 77, 032342 (2008). 

S. Maniscalco, F. Francica, R. L. ZafRno, N. Lo Gullo, 
and F. Plastina, Phys. Rev. Lett. 100, 090503 (2008). 
B. M. Garraway, Phys. Rev. A 55, 4636 (1997). 
E. Yablonovitch, Phys. Rev. Lett. 58, 2059 (1987). 
G. M. Nikolopoulos and P. Lambropoulos, J. Mod. Opt. 
49, 61 (2002). 

D. O. Guney and D. A. Meyer, J. Opt. Soc. Am. B 24, 
283 (2007). 

S. John and J. Wang, Phys. Rev. Lett. 64, 2418 (1990). 
S. John and T. Quang, Phys. Rev. A 50, 1764 (1994). 
S. M. Barnett and P. M. Radmore, Methods in Theo- 
retical Quantum Optics (Oxford University Press, New 
York, 1997). 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

R. Sprik, B. A. van Tiggelen, and A. Lagendijk, Euro- 

phys. Lett. 35, 265 (1996). 

K. Busch and S. John, Phys. Rev. E 58, 3896 (1998). 
M. Woldeyohannes and S. John, Phys. Rev. A 60, 5046 
(1999). 

E. Hagley et al., Phys. Rev. Lett. 79, 1 (1997). 

E. Ozbay et al., Appl. Phys. Lett. 64, 2059 (1994). 

P. Lodahl et al.. Nature 430, 654 (2004). 

N. H. Bonadeo et al.. Science 282, 1473 (1998). 

B. Aln, F. Bickel, K. Karraib, R. J. Warburton, and P. M. 

Petroff, Appl. Phys. Lett. 83, 2235 (2003). 

R. J. Warburton et al., Phys. Rev. B 65, 113303 (2002). 



